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ABSTRACT 



Context. The number and spatial distribution of confirmed quasi-stellar objects (QSOs) behind the Magellanic system is limited. 
This undermines their use as astrometric reference objects for studies of proper motion and of the interstellar medium along the 
line of sight. 

Aims. We search for criteria to identify candidate QSOs using near-infrared observations from the VISTA survey of the 
Magellanic Clouds system (VMC). The VMC survey provides photometry in the YJKs bands and 12 epochs in the band 
with unprecedented sensitivity and spatial resolution. 

Methods. The {Y - J) versus {J -K^) diagram has been used to distinguish QSO candidates from Milky Way stars and stars of the 
Magellanic Clouds. Then, the slope of variation in the band has been used to identify a sample of high confidence candidates. 
These criteria were developed based on the properties of 1 17 known QSOs presently observed by the VMC survey where the 
epochs span a time range up to 600 days. 

Results. YJKs magnitudes and A's light-curves of known QSOs behind the Magellanic system from present VMC data are pre- 
sented. About 75% of them show a slope of variation > lO""* mag/day and the shape of the light-curve is in general irregular and 
without any clear periodicity. A method to identify QSOs based solely on the VMC data is proposed using YJK^ colours and 
A's variability. The number of QSO candidates found in tiles including the South Ecliptic Pole and the 30 Doradus regions is 22 
and 26, respectively, with negligible contamination by young stellar objects, planetary nebulae, stars and normal galaxies. The 
high confidence in the nature of the selected objects is supported by recent studies of possible contaminants, but remains to be 
confirmed by spectroscopic follow up. 

Conclusions. By extrapolating the number of QSO candidates to the entire VMC survey area we expect to find about 1500 QSOs 
behind the LMC, 600 behind the SMC, 300 behind the Bridge and 50 behind the Stream areas. Without magnitude restrictions 
and allowing for de-blended objects, the VMC survey can find candidates > 1 mag fainter than those presented in this work. 
The VMC survey has the potential to reveal the remaining 85% of the QSOs behind the Magellanic system. Further, the A's 
light-curves can help support investigations of the mechanism responsible for the variations. The confirmed QSOs found from 
these candidates will provide a good astrometric reference grid for proper motion studies of the Magellanic system. 

Key words. Surveys - Magellanic Clouds - quasars: general - Infrared: galaxies 



1. Introduction 

The astrometric accuracy and the photometric sensitivity of 
observations made with VISTA is sufficiently good that we ex- 
pect data from the VISTA Magellanic Clouds survey (VMC; 
PI Cioni; Cioni et al. 1201 II hereafter Paper I) can be used 
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to derive proper motions of the Magellanic Clouds (MCs). 
Such proper motion studies require a reference grid of bright, 
distant, non-moving, point like objects. Quasi-stellar objects 
(QSOs) provide such a grid (e.g. Kallivayalil et al. 2006; Costa 
et al. 2011). QSOs are point-like sources believed to be pow- 
ered by accretion onto black holes in the centre of distant 
galaxies. Besides being reliable astrometric reference objects 
for studies of proper motion, they can also be used as back- 
ground sources to examine the composition of the MC inter- 
stellar medium along the line of sight (e.g. Redfield et al. 120061 
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van Loon et al. I2009I I. and are important for studies of galaxy 
formation and evolution. 

The density of QSOs with / < 19 mag is ~ 1 1 per deg^ as 
estimated from the Sloan Digital Sky Survey (SDSS) Quasar 
Catalogue (Schneider et al. 120101 ) but discovery of candi- 
date QSOs behind the MCs is complicated by the necessity 
to distinguish candidate QSOs from the dense stellar content 
of the MCs themselves. Candidate QSOs must then be ob- 
served spectroscopically to confirm which are true QSOs, and 
to make this process efficient the sample of candidate QSOs 
should be as clean as possible. Selection of candidate QSOs 
behind the MCs has been greatly improved by long-term 
multi-epoch observations, as part of micro-lensing projects 
such as the MAssive Compact Halo Objects (MACHO - 
Alcock et al. 2000) and the Optical Gravitational Lensing 
Experiment (OGLE - Udalski et al. 119921 ). and using data at 
various wavelengths, from ultraviolet (UV) to infrared (IR), 
that permit better removal of stellar objects (young stellar ob- 
jects, planetary nebulae, hot and red stars) from QSO candi- 
date samples. Methods used to identifying candidate QSOs in- 
clude X-ray (Shanks et al. 1991 ) or radio emission (Schmidt 
119681 ). mid-IR (Stern et al. 2005) and near-IR colours (see be- 
low). Also flux variations probably associated with the accre- 
tion disc have been used (Hook et al. 1994). QSO candidates 
are spectroscopically confirmed on the basis of optical and UV 
ionic emission lines, from which their redshifts are measured 
(e.g. Vanden Berk et al. 120011) . The most recent such investiga- 
tion by Kozlowski et al. ( 120121 ) focused on the southern edge 
of the Large Magellanic Cloud (LMC) and is relatively com- 
plete for objects with / < 19.2 mag, with a candidate sample 
based on Spitzer Space Telescope mid-IR colours. X-ray emis- 
sion and/or optical variability. Spectra of their sample quadru- 
pled the number of confirmed QSOs behind the LMC. 

Currently there are 360 known QSOs behind the 
Magellanic system of which 233 are behind the LMC, 100 
behind the Small Magellanic Cloud (SMC) and 27 behind the 
inter-cloud region including the Magellanic Bridge, and many 
more QSO candidates awaiting follow-up observations to es- 
tablish their nature (e.g. Kozlowski & Kochanek 2009; Kim 
et al. 120121 ). However these objects cover a limited area com- 
pared to the extent of the whole MC system being surveyed in 
the VMC survey which covers a wide area around each galaxy, 
including the Bridge region and a few fields in the Magellanic 
Stream. VMC detects sources as faint as = 23.4 mag (AB) 
with S/N= 5, coiTesponding to the luminosity of sources be- 
low the oldest turn off point of main sequence stars in the 
LMC which occurs at / ~ 22 mag. The YJK^ VMC survey, 
which is multi-epoch in K^, has the potential to considerably 
enlarge the parameter space for the search of QSOs, behind 
the Magellanic system, both in terms of sensitivity and spatial 
distribution. Near-IR criteria have been used to select high-z 
QSOs using data from the UKIRT Infrared Deep Sky Survey 
(UKIDSS) and SDSS data (Mortlock at al. 2012), and the 
highest redshift quasar currently known (z= 7.085, Mortlock 
et al. 1201 II ) was found, whilst Findlay et al. ( 120121 ) used cri- 
teria based on single epoch VISTA ZYJ data to select high 
redshift (6.5 <z< 7.5) QSOs. Kouzuma & Yamaoka (2010) 
proposed selection of candidate QSOs based on 2MASS JHK 
bands, and a series of papers have described the K excess 
(KX) method using UKI DSS J K bands and SDSS Vg bands 
(Warren, Hewett & Foltz |2000l Maddox et al. 120X21 ). 



Our aim is to establish VMC YJK^, bands selection criteria 
using the known QSOs behind the Magellanic system which 
have already been observed (by VMC), and investigate their 
utility in identifying new candidate QSOs in the areas being 
surveyed. For our purposes low redshifts are as good as high 
so the sophistication of the techniques for finding high-z QSOs 
is not needed. The first stage of the method adopted for VMC 
regions uses the (Y - J) vs {J - K^) diagram and the second 
stage uses variability. 

Section|2]describes the known QSO population behind the 
Magellanic system, and the VMC data acquired up to end 
of September 2011, before presenting the VMC photometry 
of those QSOs that have been covered so far. Section [3] uses 
VMCs YJKs photometry of these QSOs to establish colour 
criteria for isolating them from Milky Way (MW) and MC ob- 
jects, and to examine the variability of the QSOs over a 
baseline of up 12 epochs over two years. Section |4] discusses 
how the combination of colour and variability criteria perform 
on selecting out non-QSOs, using catalogues of various types 
of objects towards the MCs. Section |5] presents our conclu- 
sions on the potential of the YJKs photometry and vari- 
ability QSO selection method, and the Appendix provides the 
present sample of VMC Ks band light-curves of known QSOs. 



2. Data and Sample 

2.1. VMC data 

The near-IR data analysed in this study were obtained with the 
Visual and Infrared Survey Telescope for Astronomy (VISTA; 
Emerson & Sutherland 20101 ) for the VMC survey and in- 
clude observations acquired until the end of September 201 1. 
The data were reduced onto the VISTA photometric system 
(Vegamag = 0) with the VIST A Data Flow System (VDFS) 
pipeline v 1 . 1 (Irwin et al. 2004) and extracted from the VISTA 
Science Archival (VSA) using data release VMCv20120126. 
The VMC survey strategy involves repeated observations of 
tiles across the Magellanic system, where one tile homoge- 
neously covers an area of 1 .5 deg^ with 3 epochs in both Y 
and J filters, and 12 epochs in K^^ with the epochs spread over 
a time range of a year or longer Details about the observing 
strategy and the data reduction are given in Paper I. 

Table [1] lists the tiles which contain known QSOs and 
whose observations were by end of September 2011. Tile 
names are specified in column 1 and the limiting magnitude 
corresponding to sources with photometric errors < 0.1 mag 
in columns 2-4 while the number of available epochs is 
given in column 5 (see also Sect. 13.41 ). 

VMC sources have numerical quality flags and in the fol- 
lowing analysis we discriminates sources that are well de- 
tected (photometric flag A - with quality flags = 0-16 where 
16 indicates that a source has been de-blended in at least one 
wave band), and those with a low quality detection (photomet- 
ric flag B - with quality flags > 16 that would arise if a source 
is located in the under-exposed edge areas of a tile, or in the 
upper half of VISTA's detector #16 which is known to have a 
varying quantum efficiency, that may result in bad flat field- 
ing and unreliable magnitudes, or just a low confidence in the 
aperture magnitude of sources). See Paper I for details. 



' http ://horus .roe . ac . uk//login.html 



2 



Cioni et al.: QSOs behind the Magellanic system 



Table 1. VMC tile sensitivity and epochs. 



Limiting Magnitude" Epochs 
Tile Y J K, K, 
(mag) (mag) (mag) 

9 
10 
3 
11 
9 
1 
8 
2 
10 
10 
14 
1 
2 
12 
14 
1 



" For sources with photometric error < 0. 1 mag. 
* VMC observations are completed. 



2.2. Known quasars 

The spatial distribution of the 360 spectroscopically confirmed 
QSOs behind the Magellanic system is shown in Fig.[T] Most 
of them, confirmed prior to July 2009, are included in the 
compilation by Veron-Cetty & Veron ( 1201 Oi l and we use their 
coordinates, corrected according to Flesch ( |2012l l, instead of 
those from the original studies. The lists of QSOs behind the 
LMC and SMC are discussed in Sects. IZTTl and lzni below. 
respectively. For the presently studied VMC dataset only 117 
quasars are included within the regions observed by the VMC 
survey, with most of them contained in tiles LMC 5_5, LMC 
6.4, SMC 4_3 and SMC 5.4 (see Tab.[I]for the sensitivity lim- 
its). Tile numbering begins from the bottom right corner, in- 
creasing from right to left and from bottom to top. The first 
LMC tile is 2_3, the first SMC tile is 2_2, the first Bridge tile 
is 1^ and Stream tile 1_1 is right above the Bridge while 2_1 
is to the right of the SMC. The central coordinates of VMC 
tUes are given in Paper I while the contour plots showing the 
outer structure and bar of the galaxies as well as the distribu- 
tion of evolved giant stars are described in Cioni, Habing & 
Israel (|2000). 

Figure [T] shows that the distribution of known QSOs is in- 
homogeneous and biased to the central region of the galaxies. 
This is because most QSOs were extracted from micro-lensing 
surveys that have focused their observations on the densest 
regions, i.e. the bar of the galaxies. In addition, not all can- 
didates from these searches have been spectroscopically ob- 
served, further biasing the distribution of confirmed sources. 

2.2.1. LMC 

The list of 233 spectroscopically confirmed quasars behind 
the LMC comprises 75 QSOs from the Veron-Cetty & Veron 
(1201 Oi l compilation among which 10 were identified as X- 
ray sources with the Chandra X-ray Observatory satellite and 



3 as a result of mining OGLE-II light-curves (Dobrzycki et 
al. 120021 12005] | while 38 re fer to variability in the MACHO 
database (Geha et al. '2003'). In addition, two QSOs used by 
Anguita et al. (2000 ) and one by Pedreros et al. ( 120021 ) to study 
the proper motion of the LMC are attributed to a private com- 
munication by J. Maza in 1989 and/or do not appear in previ- 
ous studies. More recently, 1 QSO was confirmed by Hony et 
al. iiorn and 145 QSOs by Kozlowski et al. (120121 ). the lat- 
ter identified from OGLE-III light-curves. The entire MACHO 
database, with light-curves spanning a time range of ~ 7.5 yr, 
was searched by Kim et al. ( 201 II |20T2 ) who trained a sup- 
port vector machine model with diagnostic features based on 
mid-IR colours, spectral energy distribution red-shifts and X- 
ray luminosity of previously known QSOs to identify 663 high 
confidence candidates, but none of them are at present spectro- 
scopically confirmed. Finally, as part of a spectroscopic study 
on post asymptotic giant branch (AGB) stars, Kamath et al. 
(in prep.) confirmed 6 QSOs and detected one broad emission 
line in 3 additional objects that may also be QSOs. 

2.2.2. SMC 

The list of 100 spectroscopically confirmed quasars behind 
the SMC comprises 41 QSOs from the Veron-Cetty & Veron 
(|2piO) compilation among which 16 were discovered as coun- 
terparts to X-ray sources (Tinney, Da Costa & Zinnecker 
[19971 Tinney [T9991 Dobrzycki et al. 2 003bi), 3 were identi- 
fied from their OGLE-II light-curves (Ever l20()2l Dobrzycki et 
al. '2003a) and 9 from their MACHO light-curves (Geha et al. 
2003). More recently, Kozlowski et al. (201 1 ) inspected mid- 
IR colours from the Spitzer Space Telescope and the OGLE-II 
light-curves of sources in the central 1.5 deg^ of the SMC. 
They confirmed 29 QSOs and showed broad emission lines, 
typical of QSOs, in 12 additional objects but the signal to noise 
ratio (S/N) in their spectra was of lower quality. The newest 
QSOs are from Kamath et al. (in prep.) where 17 are con- 
firmed and 1 is a likely candidate since only one typical QSO 
broad emission line was detected in its spectrum. 

2.3. Cross-correlation 

VMC detections of the 126 known (spectroscopically con- 
firmed) QSOs covered by VMC data up to September 2011 
were selected by searching for the nearest counterpart within 
2". This query to the VSA produced 117 matches in the 
merged catalogue ("vmcsource") containing Y, J and band 
sources extracted from "tiledeepstacks" which are deep tile 
images resulting from the combination of individual tile im- 
ages taken at different observing epochs. A further 8 QSOs 
were matched in the un-merged catalogue, and only one 
known QSO was undetected. We inspected the images of all 
126 objects and most were good and well matched, and we 
discuss all the matches in Sect. 12.3.1! 

2.3.1. VMC-QSO distances 

Forty QSOs in the LMC from the list by Kozlowski et al. 
(2012) have counterparts at a mean distance of 0.33 + 0.17". 
Five matches were rejected because the corresponding VMC 
source, contrary to the distribution of most counterparts, was 
found at a distance > 1" and in all five cases inspection 



SMC 3.3 20.659 

SMC 3.5 21.127 

SMC 4.3 20.152 

SMC 5.4 20.678 

BRI2.3 21.222 
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LMC 4.2 21.305 

LMC 4.6 20.652 
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LMC 6.4 19.259 
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LMC 8.3 20.882 

LMC8.8'' 21.091 

LMC 9.7 20.989 



20.412 19.607 

21.002 19.693 

19.977 19.071 

20.461 19.621 

20.952 19.734 

20.742 18.418 

20.593 19.501 

20.378 18.587 

19.626 19.081 

18.962 18.576 

19.514 18.995 

20.585 17.911 

19.949 18.802 

20.364 19.676 

20.760 19.856 

20.414 18.226 
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Fig. 1. Distribution of known quasars behind the Magellanic system (top) with enlargements on the LMC (bottom-left) and SMC 
(bottom-right). All known quasars are shown with filled small black squares while those positionally matched with presently 
available VMC data are shown in red. Large empty blue rectangles indicate the VMC tiles. Contour plots show the distribution 
of evolved giant stars which delineate the outer structure and bar of the galaxies. Coordinates are given with respect to (ao, 
6oM5r, -69°). 



of the VMC images showed another, fainter, object with the 
typical red colour of QSOs present but not found in the list 
of nearby objects. A cross-correlation with the "vmcdetec- 
tion" un-merged catalogue shows that the likely counterparts 
are indeed located at distances < 1" (Table [3]). The counter- 
parts to Kamath et al.'s sources lie at a similar distance, i.e. 
0.31+0.17". 



Seventeen QSOs were matched with the sample by 
Dobrzycki et al. (120021 12003al I2003bl 120051) of which 7 in the 
SMC and 10 in the LMC. Thirty SMC QSOs were matched 
from the lists of Kozlowski et al. (201 1); 20 from their con- 
firmed list and 10 from their plausible list. Three matches were 
rejected because their detection in Y and J is suspicious; they 
are faint and close to a brighter object, and their suggested 
counterparts are at a distance > 1.3" on average. The likely 
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counterparts to these QSOs were found in the un-merged cat- 
alogue (TableO. The VMC counterpart to QSOs J0048 18.76- 
732059.6, J00483 1.50-732339.9, J010057.77-722230.8 and 
JOlOl 37.52-72041 8.9 at ~ 1" are perhaps affected by the ex- 
tended nature of the source since all other counterparts lie at 
smaller distances. Four QSOs were matched with the compila- 
tion by Veron-Cetty & Veron (12010!) of which 3 in the Bridge 
and 1 in the LMC; the latter, RXS J05466-6415, is the bright- 
est near-IR quasar in the sample with Ks - 13.590 + 0.004 
mag. These QSOs were originally confirmed by Perlman et 
al. ( fT998] l and Healey et al. ([2008b for the Bridge and LMC, 
respectively. 

The mean distance of VMC counterparts for quasars from 
Do brzyck i et al. (2002] I2003al i2003b 2005 ), Kozlowski et 
al. ( 1201 II ) and Veron-Cetty & Veron ( 120101 ) are very simi- 
lar and correspond to 0.64 + 0.30". The larger distance of 
VMC counterparts to the quasars from Kozlowski et al. ( 1201 II ) 
compared to the quasars from Kozlowski et al. (2012') is most 
likely due to their identification from the OGLE database part 
II in one case and part III in the other. OGLE-II data were 
obtained with a CCD camera with a scale of 0.4"/pixel while 
OGLE-III data were obtained with a new camera with a scale 
of 0.26"/pixel and this may introduce a systematic offset in 
the world-coordinates between the two systems. Both OGLE 
surveys are well-matched to the VISTA camera which has a 
scale of 0.34"/pixel. 

Finally, from the list of Geha et al. (|2003|) 16 QSOs, of 
which 2 in the SMC, were successfully matched. MACHO 
061.08072.0358 was excluded because two VMC candidate 
sources are present, one at 1.36" and one at 1.49". The mean 
distance of the VMC counterparts is 0.86 ± 0.29". 

With the exception of one QSO, as mentioned above, all 
known QSOs falling into VMC observed tiles were matched 
with the merged or un-merged catalogues. In the following 
analysis we consider only those found in the merged catalogue 
since the magnitudes of those in the un-merged catalogue may 
be influenced by the presence of bright neighbours. 

2.4. Summary of detected quasar sample 

Table 12] lists the information about the 117 known quasars 
found in the merged VMC catalogue. The first column gives 
the QSOs name from the literature for which references are in 
column 2. Columns 3 and 4 list the J2000 coordinates of the 
VMC counterparts. The Y, J and magnitudes and 1 sigma 
uncertainties are listed in columns 5-10, where missing values 
are indicated with -99.999. Column 1 1 lists the VDFS source 
classification flag (1 for galaxies, -3 for probable galaxies, -1 
for stars and -2 for probable stars), column 12 the VMC tiles 
where QSOs were found and column 13 the photometric flag. 

Of the 1 17 matched QSOs, 67 lie behind the LMC, 47 be- 
hind the SMC and 3 behind the Magellanic Bridge regions. 
The majority of the quasars are detected in all three wave 
bands, YJKs, with 1 QSO present only in YJ data, 3 only in 
JKs data and 2 only in A's data. The VDFS pipeline classi- 
fies ~ 53% of the matched QSOs as galaxies, among which 
two show clearly extended discs, and the other 47% as stars or 
probable stars. 

For the 8 QSOs detected in the un-merged VMC cata- 
logue, but blended with other sources which appear in the 
VMC merged catalogue. Table [3] lists the name, reference, 
wave bands and VMC tile in columns 1, 2, 3 and 4, respec- 



Table 3. VMC blended quasai-s. 



i\ alllc 


Rpf 

IvCi . 


p liter 


Ti1p 
1 lie 


J004753. 62-724350.6 


7 


K, 


SMC 4_3 


J005444.70-724813.7 


7 


K, 


SMC 4.3 


J0057 14. 13-723342.8 


7 


K, 


SMC 5.4 


MQSJ05 1944.39-701957.3 


8 


J,K, 


LMC 5.5 


MQS J052300. 14-70 183 1.7 


8 


Y,K, 


LMC 5.5 


MQS J05243 1 .88-70223 1.6 


8 




LMC 5.5 


MQS J052908.79-702445.7 


8 




LMC 5.5 


MQS J053304.31-714848.4 


8 




LMC 4.6 



References as in table[2] 



tively. These objects are listed for completeness but are not 
used in our analysis. 

3. Results 

3.1. YJKs colour-colour diagram 

3.1.1. YJKs colour-colour diagram and known QSOs 

Fig. |2] shows the distribution of the 117 known QSOs in the 
YJK^ colour-colour diagram superimposed to the distribution 
of sources detected in the LMC tile 8.8. This tile was chosen 
for comparison because it is located in the disc of the LMC 
but avoids the central region, where crowding influences the 
detection of sources rendering the sample with small photo- 
metric errors shallower than in the outer regions. In addition, 
the data for this tile are complete, i.e. all VMC epochs at each 
wave band have been obtained, reduced and included in the 
analysis. All 117 known QSOs are shown regardless of their 
quality flag and we have checked that the analysis of this paper 
does not change if only QSOs of photometric class A (Sect. |2]) 
are used. 

Most QSOs are located within the triangle limited by the 
following lines and the plotted region: 

(J-K,)^-l.25x(Y-J) + l.05 (1) 

and 

(J-K,)^2.05x(Y-J)-0.l5 (2) 
where the line 

(J-K,)^-l.25x{Y-J) + l.90 (3) 

marks the division between star-like and galaxy-like QSOs. 
This line is parallel to the first line (Eq. [U that defines a blue 
colour boundary beyond which no known QSOs are present. 
The other line (Eq.|2]) runs parallel to the locus of stars and di- 
vides them from extra-galactic sources. The selection criteria 
that embrace the sample of known QSOs are not influenced 
by the choice of tile for the foreground sources which are only 
shown for the purpose of visualisation. Note that as well as the 
QSOs there are many other objects populating this region of 
the colour-colour space (on paper the eye may be drawn to not 
notice these well due to the denser grouping outside this area 
and the colours within it). 

The highest concentrations of non-QSOs are attributed to 
main-sequence stars at the bluest colours and to giant stars of 
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Name 



Ref. 



(h 



a 
m 



s) (d 



8 

m : 



s) 



Y 

(mag) 



(Ty 

(mag) 



(mag) 



(mag) 



(mag) 



(mag) 



Class Tile Flag 



OGLE 003850.79-73 1053. 1 

J003857.50-741000.7 

J003957.65-730603.6 

J004145.04-725435.9 

J004507.50-724121.8 

J004520.77-730301.9 

J004522.89-724541.5 

J00455 1.62-725735.0 

J004556.83-725636.5 

MACHO 208. 15799. 1085 

J004736. 12-724538.2 

J004818.76-732059.6 

J004831.50-732339.9 

J004944.41-725400.5 

J004958.75-725952.2 

J005 132.67-73 14 11. 6 

J005203.27-723830.6 

J005235.03-732940.6 

J005254.15-725628.1 

J005259. 10-7245 10.9 

J005402. 17-742733.3 

J005433.79-730924.0 

J005454.05-723726.2 

J00551 1.81-723728.3 

J005517.63-712655.1 

MACHO 207.16316.446 

J005708. 12-722257.8 

OGLE 005719.84-722533.5 

J005754.26-715630.7 

J005901.04-721330.0 

J010005.10-715921.7 

J010005.70-715723.5 

J010031.06-722444.1 

J010046.02-722131.4 

J010057.77-722230.8 

J010127.75-721306.2 

OGLE010131.45-720702.6 

J010137.52-720418.9 

J010241.90-723238.7 

OGLE 010244.89-721521.7 

J010507.66-712448.4 

OGLE010522.51-715649.7 

J010650.81-713253.0 

OGLE 010656.98-722448. 1 

OGLE 010717.66-722418.3 

QJOl 16-7259 

JOl 1652.8 1-740050.8 

1WGAJ0216.6-7331 

1WGAJ0217.7-7347 

PKS0333-729 

MACHO 42.00860. 123 

MACHO 53.03360.344 

MACHO 25.03469. 117 

MACHO 25.3712.72 

J050303.97-663345.9 

MACHO 53.3970.140 

J050550.35-675017.5 

OGLE 050736.52-684751.7 

MACHO 52.4565.0356 

OGLE 050833.29-685427.5 



3 

10 
7 
7 
7 
7 
7 
7 
7 
2 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 

10 

7 

7 

7 

10 

2 

7 

4 

10 

7 

7 

7 

7 

7 

7 

7 

4 

7 

7 

3 

10 
4 
10 
4 

4 

9 

10 

6 

6 

6 

2 

2 

2 

2 

10 

2 

5 

1 

2 

1 



00:38:50.80 
00:38:57.53 
00:39:57.64 
00:41:45.04 
00:45:07.55 
00:45:20.84 
00:45:22.95 
00:45:51.64 
00:45:56.80 
00:47:15.92 
00:47:36.13 
00:48:18.74 
00:48:31.50 
00:49:44.42 
00:49:58.86 
00:51:32.66 
00:52:03.30 
00:52:34.94 
00:52:54.23 
00:52:59.08 
00:54:02.24 
00:54:33.82 
00:54:53.97 
00:55:11.78 
00:55:17.72 
00:55:34.67 
00:57:08.13 
00:57:19.82 
00:57:54.33 
00:59:00.94 
01:00:05.13 
01:00:05.70 
01:00:31.05 
01:00:46.02 
01:00:57.77 
01:01:27.81 
01:01:31.45 
01:01:37.50 
01:02:41.90 
01:02:44.90 
01:05:07.68 
01:05:22.54 
01:06:50.90 
01:06:56.99 
01:07:17.73 
01:16:33.51 
01:16:52.85 
02:16:44.57 
02:17:45.11 
03:32:42.95 
04:46:10.96 
05:00:53.98 
05:01:46.72 
05:02:53.73 
05:03:04.0 
05:04:35.99 
05:05:50.40 
05:07:36.44 
05:08:30.62 
05:08:33.21 



-73:10:53.4 
-74:10:00.9 
-73:06:03.6 
-72:54:35.6 
-72:41:21.5 
-73:03:01.6 
-72:45:41.4 
-72:57:35.0 
-72:56:36.3 
-72:41:12.0 
-72:45:38.1 
-73:21:40.1 
-73:23:40.1 
-72:54:00.3 
-72:59:51.9 
-73:14:11.7 
-72:38:30.6 
-73:29:40.9 
-72:56:28.1 
-72:45:10.8 
-74:27:33.4 
-73:09:24.1 
-72:37:26.5 
-72:37:28.4 
-71:26:55.2 
-72:28:34.3 
-72:22:57.8 
-72:25:34.0 
-71:56:30.8 
-72:13:30.0 
-71:59:21.5 
-71:57:23.4 
-72:24:44.4 
-72:21:31.2 
-72:22:31.0 
-72:13:06.4 
-72:07:03.1 
-72:04:19.1 
-72:32:38.8 
-72:15:22.0 
-71:24:48.4 
-71:56:49.9 
-71:32:53.0 
-72:24:48.4 
-72:24:18.3 
-72:59:49.7 
-74:00:50.8 
-73:31:40.5 
-73:47:23.3 
-72:49:04.4 
-72:05:10.0 
-66:44:00.0 
-67:32:40.1 
-67:25:45.3 
-66:33:46.0 
-66:24:15.9 
-67:50:18.0 
-68:47:51.6 
-67:02:28.9 
-68:54:27.3 



16.562 
17.371 
18.321 
17.662 
19.482 

-99.999 
19.934 
19.169 
17.617 
18.685 
18.864 
19.419 
19.419 

-99.999 
19.528 
19.078 
19.948 
20.235 
19.257 
19.812 
18.243 
20.171 
19.402 
19.241 
18.325 
17.335 
19.398 
17.151 
18.133 
19.486 
19.108 
18.900 

-99.999 

-99.999 
19.492 
19.437 
18.308 
18.175 
17.885 
17.925 
17.933 
18.334 
17.871 
18.489 
19.674 
17.719 
17.455 
18.212 
19.038 
18.648 
16.522 
18.563 
17.122 
17.332 
15.395 
17.065 
16.995 
18.204 
17.856 
17.977 



0.006 
0.010 
0.022 
0.013 
0.057 

-99.999 
0.085 
0.044 
0.013 
0.029 
0.034 
0.053 
0.053 

-99.999 
0.060 
0.041 
0.088 
0.112 
0.047 
0.080 
0.017 
0.104 
0.054 
0.047 
0.017 
0.009 
0.036 
0.008 
0.015 
0.039 
0.029 
0.024 

-99.999 

-99.999 
0.039 
0.037 
0.017 
0.015 
0.013 
0.013 
0.013 
0.017 
0.013 
0.019 
0.045 
0.010 
0.010 
0.015 
0.024 
0.020 
0.005 
0.018 
0.009 
0.010 
0.003 
0.007 
0.008 
0.040 
0.014 
0.033 



16.142 
17.104 
18.334 
16.859 
19.239 
19.968 
19.767 
18.781 
17.210 
18.358 
18.514 
18.910 
18.910 
-99.999 
19.196 
18.673 
19.451 
19.331 
18.546 
19.520 
17.726 
19.614 
18.849 
18.454 
17.741 
17.454 
18.608 
16.463 
17.435 
18.915 
18.870 
17.816 
17.890 
18.600 
19.005 
19.025 
17.964 
17.463 
17.875 
17.707 
17.347 
17.838 
17.749 
17.922 
19.443 
17.155 
16.964 
17.995 
18.655 
18.496 
15.812 
18.178 
16.505 
16.899 
14.806 
16.749 
16.272 
18.168 
17.564 
17.816 



0.006 
0.010 
0.026 
0.009 
0.054 
0.100 
0.085 
0.037 
0.011 
0.026 
0.030 
0.041 
0.041 

-99.999 
0.052 
0.034 
0.066 
0.059 
0.031 
0.071 

0.0148 
0.075 
0.040 
0.029 
0.014 
0.012 
0.025 
0.007 
0.011 
0.030 
0.030 
0.015 
0.015 
0.024 
0.033 
0.033 
0.016 
0.012 
0.015 
0.013 
0.011 
0.015 
0.014 
0.016 
0.045 
0.009 
0.008 
0.015 
0.021 
0.024 
0.005 
0.021 
0.009 
0.011 
0.003 
0.009 
0.008 
0.052 
0.018 
0.038 



14.878 


0.006 


-\ 


SMr4 3 




1(S 44Q 


0193 


_1 








OIQ 


1 


CMP 4 ^ 




1 5 nsn 


006 


1 


sivrr4 3 




18.157 


0.048 


_2 


SMC 4_3 




18.709 


0.074 


1 


SMC 4 3 


A 


18.961 


0.093 


1 


SMC 4 3 


A 


17.345 


0.025 


-I 


SMr4 3 


A 


16.600 


0.015 


1 


SMC 4_3 


A 


17.747 


035 


_1 


SMr 4 3 


A 


16.842 


010 


1 


SMr 4 3 


A 


17 439 


097 


_1 




A 
/ \ 


17 439 


097 


_1 


SMr4 3 


A 


17.201 


0.023 


_1 


SMC 4 3 


A 


1 8 175 


048 


1 


<?Mr4 3 


A 


17 3Q6 


097 


_1 


SMr 4 3 


A 


17.898 


0.039 


_1 


SMr 4 3 


A 


18 518 

1 O.J 1 o 


0(S4 


_1 


SMr 4 3 


g 


17.135 


0.022 


1 


SMr4 3 


A 




OS^ 


1 
i 


SMP4 ^ 


A 


16 14Q 


01 1 


1 


SMr 3 3 


g 


18.144 


048 


_1 


SMr4 3 


A 


18.195 


0.050 


1 


SMC 4_3 


A 


17.040 


0.021 


_1 


SMC 4 3 


A 


16.805 


0.015 


_1 


SMr5 4 


g 


1(S n55 


010 


_1 


SMr 5 4 


A 


17 088 


018 


1 


SMr 5 4 


A 


14.728 


0.006 


1 


SMr 5 4 


A 


15.531 


0.008 


I 


SMr5 4 


A 


17.715 


0.026 


\ 


SMr 5 4 


A 


1 8 309 


038 


1 




A 


16 933 


016 


_1 




A 

/ \ 


17 066 


018 


_2 




A 


17.134 


0.018 


1 




A 


17.445 


0.022 


1 




A 


1 8 040 


0.032 


1 




A 


17.476 


0.022 


_1 


SMr 5 4 


A 


15.783 


000 


1 


SMr 5 4 


A 


16.278 


019 


_1 


SMr 5 4 


g 


1(S 589 


014 


_1 


SMr 5 4 


A 


15 710 


000 


1 


SMr5 4 


g 


16 803 


016 


_1 


SMr 5 4 


A 
/ \ 


16 540 


013 


_1 


SMr5 4 


A 


16.299 


0.012 


1 


SMC 5 4 


A 


18.240 


0.037 


I 


SMr 5 4 


A 


1 5 603 


01 1 


I 


SMr 4 5 


A 


16 0Q3 


0.111 


_1 


SMr 3 5 


A 


17.155 


0.019 


_1 


BRI2 3 


A 


17.153 


0.019 


I 


BRI2 3 


A 


1 7 40Q 


048 


_i 


DDT Q 7 


A 
/ \ 


14.952 


0.006 


-1 


LMC4^ 


A 


17.184 


0.019 


-1 


LMC 8_3 


B 


14.561 


0.005 


-1 


LMC7.3 


A 


15.300 


0.008 


-1 


LMC 7_3 


A 


13.173 


0.003 


1 


LMC 8_3 


A 


15.696 


0.009 


-1 


LMC 8_3 


A 


14.772 


0.006 


1 


LMC 7_3 


A 


16.729 


0.023 


-1 


LMC 6_4 


B 


16.268 


0.015 


-1 


LMC 7.3 


A 


16.345 


0.017 


1 


LMC 6.4 


A 



6 



Table 2. continued. 



Cioni et al.: QSOs behind the Magellanic system 



Name 



Ref. 



(h 



a 
m 



s) 



6 

(d : m : 



s) 



Y 

(mag) 



o-y 
(mag) 



J 

(mag) 



o-j 
(mag) 



(mag) 



(mag) 



Class Tile Flag 



MQS J050842.68-690959.4 
MACHO 20.4678.600 
J050900.76-672114.1 
J050915.35-695417.3 
OGLE 050924.05-672 1 24. 1 
J050940.60-685320.6 
MACHO 5.4892. 1971 
MQS J05 1 045 . 85-694 1 26. 5 
MQS J051049.16-692851.9 
MQS J05 1159.38-694726.1 
J05 1244.64-693857.2 
MQS J05 1345.50-690636.9 
MQS J05 1 404. 8 1 -694446. 6 
J05 1433.58-693245.7 
MQS J05 1534.35-685422.2 
MQS J05 1547.78-700002.3 
MACHO 78.5855.788 
MACHO 2.5873.82 
MQS J05 1 7 1 6.95-704402.0 
MQS J05 1739.83-695705.3 
MQS J05 1746.88-685458.1 
MQS J051834.75-703751.3 
OGLE 051853. 19-690217.7 
MQS J051856.04-701743.6 
MQS J051912.33-710914.6 
MQS J05 1 9 1 6. 98-70445 8. 8 
MQS J05 1953.65-704622.7 
MQS J052037.03-705 100.1 
MACHO 6.6572.268 
J052141.75-703029.2 
MACHO 13.6805.324 
MACHO 13.6808.521 
MQSJ052312.95-701529.8 
MQS J052335. 33-710758.1 
MQSJ052336.58-710033.9 
MQS J052354.44-704821.9 
MQS J052408. 88-704812.8 
MQS J052444.07-702601.9 
MQS J052502.45-705023.0 
MQS J052521.42-710932.2 
MQS J052528.91-700448.6 
MQSJ052529.65-701018.9 
MQS J052705. 98-702610.7 
MQS J052745.27-700533.9 
MQS J052747.29-71 1058.9 
MQS J052749.08-703641.7 
MQS J052931. 17-702958.7 
MQSJ053044.94-712358.5 
MQS J053120.33-715326.0 
MQS J053225.83-713328.6 
MQS J053708.16-721747.5 
MQS J053741. 18-715945.8 
MQS J054200.63-71401 1.4 
MQS J054227.27-720454.7 
MQS J054319.60-715151.8 
RXS J05466-6415 
MACHO 75. 13376.66 



2 
5 

10 
1 
5 
2 



05:08:42.67 
05:08:54.13 
05:09:00.67 
05:09:15.37 
05:09:23.77 
05:09:40.53 
05:10:32.39 
05:10:45.84 
05:10:49.15 
05:11:59.33 
05:12:44.54 
05:13:45.45 
05:14:04.80 
05:14:33.45 
05:15:34.35 
05:15:47.77 
05:16:26.31 
05:16:28.84 
05:17:16.93 
05:17:39.80 
05:17:46.88 
05:18:34.74 
05:18:53.08 
05:18:56.05 
05:19:12.29 
05:19:16.96 
05:19:53.63 
05:20:37.01 
05:20:56.97 
05:21:41.61 
05:22:47.18 
05:22:47.68 
05:23:12.94 
05:23:35.33 
05:23:36.58 
05:23:54.42 
05:24:08.87 
05:24:44.07 
05:25:02.47 
05:25:21.42 
05:25:28.82 
05:25:29.65 
05:27:05.89 
05:27:45.26 
05:27:47.29 
05:27:49.09 
05:29:31.19 
05:30:44.92 
05:31:20.33 
05:32:25.74 
05:37:08.15 
05:37:41.18 
05:42:00.63 
05:42:27.33 
05:43:19.59 
05:46:41.84 
06:02:34.32 



-69:09:59.5 
-67:37:34.2 
-67:21:14.1 
-69:54:17.3 
-67:21:24.5 
-68:53:20.3 
-69:27:15.7 
-69:41:26.4 
-69:28:51.8 
-69:47:25.8 
-69:38:57.1 
-69:06:37.1 
-69:44:46.5 
-69:32:45.6 
-68:54:22.1 
-70:00:02.2 
-69:48:19.2 
-68:37:02.3 
-70:44:02.2 
-69:57:05.2 
-68:54:58.0 
-70:37:51.2 
-69:02:17.5 
-70:17:43.7 
-71:09:14.6 
-70:44:58.8 
-70:46:22.6 
-70:51:00.1 
-70:24:52.7 
-70:30:28.8 
-71:01:30.7 
-70:47:34.8 
-70:15:29.7 
-71:07:58.1 
-71:00:33.9 
-70:48:21.9 
-70:48:12.8 
-70:26:01.9 
-70:50:23.1 
-71:09:32.2 
-70:04:48.3 
-70:10:18.9 
-70:26:11.2 
-70:05:33.9 
-71:10:59.0 
-70:36:42.1 
-70:29:58.7 
-71:23:58.5 
-71:53:26.0 
-71:33:28.5 
-72:17:47.5 
-71:59:45.7 
-71:40:11.4 
-72:04:54.4 
-71:51:51.8 
-64:15:22.0 
-68:30:41.5 



18.723 
18.974 
18.704 
16.721 
17.793 
17.934 
18.106 
18.097 
19.397 
18.894 
17.605 
18.615 
18.690 
17.466 
19.054 
18.798 
17.360 
16.052 
17.332 
18.883 
18.543 
19.152 
19.367 
18.670 
18.850 
18.187 
19.323 
18.737 
17.425 
18.281 
18.160 
18.175 
18.037 
19.257 
18.910 
-99.999 
18.834 
18.757 
18.967 
19.292 
17.434 
18.918 
18.522 
17.973 
19.267 
17.843 
19.048 
18.527 
18.054 
19.153 
17.542 
18.546 
18.741 
19.319 
18.874 
15.522 
18.127 



0.064 
0.032 
0.027 
0.012 
0.014 
0.032 
0.036 
0.036 
0.114 
0.074 
0.024 
0.057 
0.061 
0.021 
0.086 
0.068 
0.019 
0.007 
0.013 
0.073 
0.055 
0.055 
0.113 
0.036 
0.042 
0.024 
0.063 
0.037 
0.013 
0.026 
0.023 
0.024 
0.021 
0.057 
0.043 
-99.999 
0.041 
0.039 
0.045 
0.060 
0.014 
0.044 
0.032 
0.020 
0.059 
0.018 
0.049 
0.032 
0.015 
0.032 
0.011 
0.021 
0.024 
0.036 
0.026 
0.004 
0.016 



18.423 
18.642 
18.144 
16.054 
17.110 
17.465 
17.652 
17.676 
18.677 
18.257 
16.855 
17.866 
17.998 
17.100 
18.221 
18.525 
17.135 
15.511 
16.727 
18.623 
18.420 
18.300 
18.757 
17.588 
18.185 
17.726 
19.095 
18.111 
17.031 
18.137 
17.731 
17.708 
17.418 
18.293 
18.234 
-99.999 
18.325 
18.234 
18.705 
18.885 
17.067 
18.602 
18.076 
17.765 
18.487 
17.311 
18.549 
18.175 
17.714 
18.921 
17.138 
18.139 
18.554 
18.570 
18.658 
15.147 
17.924 



0.064 
0.035 
0.025 
0.009 
0.013 
0.027 
0.032 
0.033 
0.079 
0.054 
0.016 
0.039 
0.043 
0.020 
0.053 
0.069 
0.021 
0.006 
0.010 
0.076 
0.065 
0.033 
0.085 
0.019 
0.030 
0.021 
0.064 
0.028 
0.013 
0.029 
0.020 
0.021 
0.017 
0.032 
0.031 
-99.999 
0.034 
0.031 
0.045 
0.053 
0.013 
0.043 
0.028 
0.022 
0.038 
0.015 
0.041 
0.030 
0.014 
0.035 
0.010 
0.019 
0.025 
0.025 
0.027 
0.004 
0.017 



17.441 
17.403 
17.130 
4.482 
15.636 
16.760 
16.736 
16.381 
16.962 
16.912 
15.143 
16.750 
16.280 
16.058 
16.454 
17.628 
15.860 
14.247 
15.348 
17.285 
16.921 
17.018 
17.631 
16.775 
-99.999 
16.894 
18.041 
16.494 
16.305 
17.015 
16.894 
16.819 
15.919 
17.059 
17.401 
17.497 
17.291 
17.222 
17.658 
17.761 
16.252 
17.746 
17.155 
16.662 
16.742 
16.202 
17.145 
16.978 
16.744 
17.483 
15.502 
16.969 
17.448 
16.880 
17.120 
13.590 
16.615 



0.039 
0.033 
0.027 
0.005 
0.010 
0.023 
0.022 
0.017 
0.026 
0.026 
0.008 
0.023 
0.016 
0.013 
0.018 
0.046 
0.012 
0.005 
0.008 
0.035 
0.026 
0.021 
0.046 
0.018 
-99.999 
0.019 
0.043 
0.015 
0.013 
0.021 
0.019 
0.018 
0.010 
0.021 
0.027 
0.029 
0.025 
0.024 
0.032 
0.035 
0.013 
0.035 
0.023 
0.017 
0.017 
0.012 
0.023 
0.021 
0.025 
0.044 
0.010 
0.029 
0.041 
0.026 
0.031 
0.004 
0.034 



LMC 6. 
LMC 7. 
LMC 7. 
LMC 6. 
LMC 7. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 6. 
LMC 5. 
LMC 6. 
LMC 6. 
LMC 5. 
LMC 6. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 5. 
LMC 4. 
LMC 4. 
LMC 4. 
LMC 4. 
LMC 4. 
LMC 4. 
LMC 4. 
LMC 9. 
LMC 6. 



A 
B 
B 
B 
B 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
B 
A 
A 
A 
B 
A 
B 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
B 
B 
B 
A 
A 
A 
A 
A 
A 
A 
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Fig. 2. Colour-colour diagram showing the known QSOs 
where in red (filled squares) are those classified as galaxy-like 
and in blue (filled triangles) those classified as star-like by the 
VDFS pipeline. All sources from the LMC tile 8_8 with pho- 
tometric eiTors < 0. 1 mag and quality flags = in each wave 
band are shown with black dots. Dashed black lines identify 
the regions where known QSOs are found. In particular, re- 
gion A is dominated by QSOs with a stellar morphology and 
region B by QSOs with a detectable host galaxy, i.e. galaxy 
morphology. The LMC and MW stellar populations dominate 
region D. The green continuous line encloses the region iden- 
tified in Paper II to contain PNe. 



the LMC. The plume of sources with (Y-J) > 0.5 mag be- 
longs to the MW foreground. Source at (Y-J) ~ 0.6 mag and 
(J - Kfi) ~ 1.5 mag are mostly background galaxies. It is inter- 
esting to note that QSOs morphologically classified as galaxy- 
like from VMC data dominate region B while those classified 
as star-like are mostly located in region A blueward from it. In 
addition, many star-like QSOs are redder than the green sloped 
line delimiting the region where most planetary nebulae (PNe) 
are found (Miszalski et al, 2011a hereafter Paper II). Only a 
handful of QSOs occupy region C. The number of sources in 
this region does increase in the presence of interstellar redden- 
ing like in the case of tile LMC 6_6 (Table |4]l. QSOs located 
outside regions ABC, i.e. in region D, cannot be distinguished 
from LMC and MW stars by their YJK^^ colours. Their num- 
ber is, however, small. It would be worthwhile to explore the 
presence of QSOs in the region resulting from the intersection 
of region D and that defined by PNe. 

Based on the number of known QSOs distributed in re- 
gions A&B (Fig. |2]l, 65% and 21% of objects classified as 
star-like occupy these regions respectively, whilst the relative 
percentages of galaxy-like QSOs are 35% and 79%. Outside 
regions A&B lie ~ 16% of the QSOs, with half classified as 
star-like and half as galaxy-like. 

The YJKs colour-colour diagram (Fig. |2]i represents 
the best VMC diagram to distinguish the bulk of back- 




Mrk231 



C 0.5 1 

Y-J 

Fig. 3. Colour-colour diagram showing the distribution of 
sources in tile LMC 8_8 with photometric errors < 0.1 
mag and quality flags = in each wave band (black dots). 
Superimposed are models from the SWIRE template library 
(Polletta et al. 2007). The green and dashed black lines are as 
in Fig.|2l Most QSOs lie in regions A and B. 



ground QSOs from the foreground stellar population of the 
Magellanic Clouds and the MW. It also provides a clear indi- 
cation of the location of normal galaxies that contaminate the 
QSO sample. We note that the first part of the KX method to 
select QSOs (Maddox et al. |2012i ) instead involves gJK, and 
in a sense our step is similar except using Y instead of g, not 
least because we have Y, but not g, data of appropriate depth. 
For a similar distribution in the VJK colour-colour diagram 
see Warren, Hewett & Foltz ( l2000b . 

3.1.2. YJKs colour-colour diagram and QSO model 
templates 

QSO templates obtained from Spitzer-space-telescope Wide- 
field InfraRed Extragalactic (SWIRE) project template library 
by Polletta et al. ( i2007j have been convolved with the VISTA 
spectral response over each of the YJKf^ filters and the location 
of these QSOs is shown in Fig. [3] superimposed on the distri- 
bution of sources in tile LMC 8_8 as in Fig. |2] The average 
redshift from the QSO models (excluding 119254 and Sey 1.8) 
is 1.22 + 0.25 in region A and 0.44 + 0.25 in region B. This is 
consistent with the trends presented in Sect. 13. 31 

The template models are semi-empirical, i.e. they have 
been constructed combining observational data for similar ob- 
jects. Those shown in Fig. [3] refer to three type 1 QSOs, one 
type 2 QSO and two obscured QSOs as well as a moderately 
luminous AGN. The QSOl templates were derived from the 
combination of SDSS spectra and rest-frame IR data where 
QSOl, TQSOl and BQSOl refer to the average flux of all, 
the highest 25% and the lowest 25% of the measurements, re- 
spectively. The QS02 template was obtained by combining 
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a optical/near-IR spectrum with the rest-frame IR data. The 
AGN template, Seyfert 1.8, is the result of broad-band photo- 
metric data and spectra from the Infrared Space Observatory 
while the templates Mrk 231 and IRAS 19254-7245 South 
(119254) are fits to their composite, AGN plus starburst, spec- 
tral energy distributions. 

Regions A+B where most known QSOs are found (Fig.|2]i 
are well matched by the models and a few QSOs may be found 
also in region C. 

3.1.3. Correction for reddening 

The QSOs and stars in the Magellanic Clouds as well as stars 
in the MW are reddened by interstellar dust present along the 
line of sight and the effect of reddening has been checked as 
follows. Reddening values in terms of E{V - 1) were extracted 
from the extinction map from Haschke et al. ( 1201 11 ) for all 
but 10 (located outside the map) of the 117 QSOs. Using the 
Cardelli, Clayton and Mathis (1989) extinction law, the aver- 
age QSO extinction is A{Y) = 0.045 mag, A{J) = 0.032 and 
A(Ks,) = 0.013 mag. The estimated reddening, while account- 
ing for the MW and LMC components, does not account for 
the intergalactic dust that lies between the LMC and the QSOs 
and for dust in the torus surrounding the AGN. 

The stellar population shown in Fig. |2] corresponds to the 
LMC and MW stars located in tile LMC 8_8. This tile is out- 
side the region studied by Haschke et al. $201 IK Therefore, 
reddening values were extracted from the Zaritsky et al. ([2004) 
extinction map. Here, instead of de-reddening each individual 
star, an average value referring to cool stars only, A(V) - 0.35 
mag, was obtained in a region centred on the tile and with a 
radius of 12'. Cool stars were chosen because they have simi- 
lar characteristics to red clump giant and RR Lyrae stars used 
by Haschke et al. (1201 11 1. The size of the region, whilst be- 
ing smaller than the extent of a tile, is representative of the 
average extinction since this tile is located in the outer disc 
of the LMC where the reddening is uniform. By applying the 
extinction law from Cardelli et al. (119891) extinction values of 
A(F) = 0.14 mag, A(J) = 0.10 mag and A(K,) = 0.04 mag 
were derived. These values may be too high since Haschke et 
al. ( 1201 l i ) found that in the less populated regions of the LMC 
differences in E{V - I) between their extinction map and that 
of Zaritsky et al. ("2004) amount to 0.07 - 0.10 mag on av- 
erage. In addition, while these values are applicable to LMC 
stars they are certainly too high towards foreground MW stars 
and too low towards extra-galactic sources. 

In view of the uncertainty in the extinction of LMC tile 8_8 
outlined in the previous paragraph, the Schlegel, Finkbeiner & 
Davis ( 1998) extinction map based on DIRBE/IRAS mid-IR 
observations was used to find average reddening and extinc- 
tion, E(B -V) = 0.067 and A(V) = 0.207 mag, in a region of 
2 deg- in size centred on tile LMC 8_8. The extinction value is 
consistent with the mean value, A(y) = 0.187+0. 040 mag, de- 
rived in Paper IV from the recovering of the star formation his- 
tory. A(Y) = 0.083 mag, A(7) = 0.058 mag and A(K,) = 0.024 
mag are then obtained. These values are 0.02 - 0.06 mag 
smaller than those derived from the Zaritsky et al. (2004) ex- 
tinction map. Using the TRILEGAL code (Girardi et al. 2005 ) 
to simulate the MW population we obtain that it accounts for 
A(V) - 0.05 of absorption implying that most of the reddening 
towards the LMC is due to the galaxy itself. 



The un-reddened colours of QSOs and stars show that the 
boundaries are not significantly affected and the standard devi- 
ation of the reddening values was about as high as the values 
themselves, so we do not apply any reddening correction in 
our analysis. 

3.2. Colour-magnitude diagrams 

The distribution of known QSOs in the VMC colour- 
magnitude diagrams (CMDs) is shown in Fig.|4]superimposed 
on the distribution of LMC and MW sources present in tile 
LMC 8_8. 

The (Y - J) vs. F or y CMDs show that QSOs overlap 
with the stellar populations of the LMC and of the MW that 
is also not distinguishable from LMC stars. The {Y - Ks) vs. 
Y or CMDs provide the best separation between the LMC 
and MW stellar population as well as background galaxies. 
QSOs, however, especially those with a stellar morphology, do 
overlap with MW stai's. The (J - K^) vs. J or CMDs offer 
instead the best separation between QSOs and stars of the MW 
or LMC. In the (J - K,) vs. CMD, galaxy-like QSOs do 
overlap with the cone of galaxies (with base 1 < (7 - K^) < 2 
mag and vertex at (/-/Ts) ~ 1 .5 mag); see Kerber et al. ( 120091 ), 
while star-like QSOs seem well separated from it. 

The distribution of known QSOs in the CMDs unsurpris- 
ingly indicates that the present sample is magnitude limited 
and that the VMC data potentially improves the magnitude 
limit for finding candidates fainter than the present known 
sample QSOs by > 1 mag. 

3.3. Trends witli redshift 

Figure |5] shows the distribution of redshifts with respect to 
VMC magnitudes and colours. There appears to be no trend 
as a function of magnitude nor of (Y - J) colour The QSOs 
occupy the whole of the parameter space delimited on one side 
by the bright limit of the survey and on the other by the faint 
limit of the present sample. The faint limit is also influenced 
by the incomplete observations of epochs. 

The plots of redshift as a function of ( J - /Tj) and (Y - Ks) 
colours suggest a trend of larger redshifts at smaller colours. 
This indicates that QSOs at small redshifts are redder than 
QSOs at larger redshifts as predicted by Maddox & Hewett 
( 120061 ). These authors simulated the evolution of quasars in the 
near-IR and found them getting bluer as redshift goes from 
to 3. Figure |5] also shows a dichotomy between star-like and 
galaxy-like QSOs where the latter are concentrated at small 
redshifts as in Maddox et al. (120081 ). 

3.4. Variability from photometry 

The VMC survey provides at least 12 epochs in observed 
under similar observing conditions (Paper I). Additional 
epochs may, however, be obtained in worse conditions. The 
number of epochs currently available (Table[Tli is not sufficient 
to sample in detail a QSO light-curve but the accuracy of the 
photometry is enough to detect significant variability. 

The light-curves for each known QSO from the 
presently available VMC data are shown in the Appendix. 
Measurements obtained during the same night, excluding 
those that do not meet the observing constraints have been av- 
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Fig. 4. Hess diagrams of sources in tile LMC 8_8 with photometric errors < 0.1 mag and quality flags = 0, in each pair of two 
wave bands, shown as contour plots. The known QSOs behind the Magellanic system are displayed as squares. QSOs are colour 
coded as in Fig. |2] 



eraged and the error bar corresponds to the standard deviation 
of the mean. We refer to the ~ 39% / ~ 27% / ~ 33% of 
the known QSO light-curves sampled over 300 - 600/40 - 
80/shorter days as the long /mid /short-range samples. The 
current short sample includes some single epoch observations. 
The acquisition of band epochs for most VMC tiles is still 
on-going (Table [1). 

An example of QSO /Ts light-curve is shown in Fig|6] The 
shapes of the light-curves are in general irregular without any 
clear periodicity. There is often a changing brightness as a 
function of time both in the long- and mid-range samples. 
Peaks superimposed on this trend are in some cases significant 
with respect to the error bars. On average there are broad (up to 
100 days wide) and narrow (20 - 40 days wide) peaks. Curves 
of the mid-range sample appear noisier but this is largely an 
efi'ect of the more frequent sampling compared to the long- 
range sample. Long-range curves do show a similar behaviour 
in the time range where more data points are available. 

To quantify the variation the slope in the Kf^ band has been 
derived for each known QSO and the number distribution of 
the slope values is shown in Fig.|7] There is a narrow peak due 
to sources with Slope ~ mag/day while most of the sample 
has 1 < Slope < 3.5 x 10""^ mag/day. Only with a larger sample 
of QSOs we will be able to establish if the latter is consistent 
with a decreasing gaussian function or represents a bump in 
the distribution. The remaining QSOs define a constant tail out 



to Slope ~ 10" mag/day. In summary ~ 75% of the presently 
known QSOs have Slope > 10 '* mag/day. In the next section 
we will explore the potential of using variability, along with 
the colour regions given above, as a discriminant for candidate 
QSOs. 

3.5. Selection of candidate QSOs 
3.5.1 . Expected number of QSOs 

In a given 1.5 deg^ V MC tile we expect 17 QSOs with / < 19 
mag (Schneider et al. 120101 ). The approximate ; band mag- 
nitude corresponding to a given Ks band magnitude has been 
calculated assuming that the ~ 2.5 mag difference between the 
tip of the red giant branch (TRGB) at and at / is applica- 
ble to other sources and that there is not a major difference 
between the / and / filters. The TRGB in the LMC occurs at 

= 12 mag and at / = 14.54 mag (Cioni et al. 2000). The 
VMC survey, however, can probe fainter magnitudes (see Fig. 
|4]) and the number of QSO candidates is therefore larger 

Tile LMC 5_5 is entirely within the area where the known 
QSO sample is complete to / < 19.2 mag (Kozlowski et 
al. 120121 ) and in this tile the faintest known QSO (MQS 
J05 1953.65 - 704622.7) has VMC magnitudes Y = 19.32, 
J = 19.09, Ks = 18.04 mag. The numbers of known QSOs 
in regions A, B and C of tile LMC 5_5 are listed in Table |4] 
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Fig. 5. Distribution of redshift versus magnitudes (top line) and colours (bottom line) for known QSOs behind the Magellanic 
system. Filled symbols indicate sources classified as star-like by the VDFS pipeline and empty symbols those classified as 
galaxies. Note the possible trend of redshift as a function of (7 - A^s) and (Y - K^) colours, and the concentration of galaxy type 
objects with small redshift at red colours. 



along with related information for the two other completed 
tiles 6_6 and 8_8, together with the number of objects brighter 
than these limits in each tile and with photometric errors < 0.1 
mag and quality flags = in each wave band. Even with the 
small number statistics it is clear that galaxy-like QSOs dom- 
inate region B and star-like QSOs dominate region A in these 
three tiles. 

By taking the ratio of all sources lying in colour selected 
regions A & B (Fig.|2l) found down to the same magnitude as 
the faintest and with the same photometric and quality flags 
to the number of known QSOs, it appears that ~ 6% of all 
sources in tile LMC 5_5 are known QSOs. Those in region 
C, may be influenced by reddened LMC sources and so their 
percentage should be considered as an upper limit. Using this 
information it is possible to extrapolate the number of QSOs 
VMC may expect to find in the other tiles. 

Using Table HI we now analyse the other two LMC tiles 
for which VMC observations are completed. Tile LMC 6_6 
includes the 30 Doradus regions where the reddening is the 
highest within the LMC implying the presence of many red- 
dened LMC sources in region C. Tile LMC 6_6 contains more 
sources than tile LMC 8_8, on average, because of it is nearer 
to the centre of the galaxy, but the difficulty in detecting them 
clearly results in larger photometric errors and quality flags. 
This explains the inferior number of sources selected to have 
photometric errors < 0.1 mag and quality flags = in each 
wave band. Taking all objects in colour regions A & B and as- 
suming the same percentage of known QSOs, as derived from 
tile LMC 5_5, we find the number of QSOs VMC might expect 
to detect is 13 in tile LMC 6_6 and 30 in tile LMC 8_8. 



Table 4. VMC quasar candidates. 



Sources 


Tile 


Region A" 


Regions" 


Region C° 






gal. + star 


gal. + star 


gal. + star 


No. QSOs 


LMC 5.5 


4-H8 


5 + 


2+1 


Mag. Lim. 


LMC 5.5 


137 + 29 


140 + 3 


14 + 21 


QSOs % 


LMC 5.5 


2.9 + 27.6 


3.6 + 


14.3 + 4.8 


Mag. Lim. 


LMC 6.6 


86 + 28 


97 + 7 


152* + 845* 


Exp. No. QSOs 


LMC 6.6 


2 + 8 


3 + 


22* + 40* 


Var. QSOcand. 


LMC 6.6 


15 


11 


49* 


Mag. Lim. 


LMC 8.8 


129 + 55 


317 + 1 


3 + 12 


Exp. No. QSOs 


LMC 8.8 


4+15 


11+0 


0+ 1 


Var. QSOcand. 


LMC 8.8 


14 


7 


1 



" Regions ABC are those in Fig.O 
* Includes reddened LMC stars. 

Mag. Lim. corresponds to the No. objects with Y < 19.32 mag, 
J < 19.09 mag and < 18.04 mag, photometric errors < 0.1 mag 
and quality flags = in each wave band. 



3.5.2. Candidate QSOs from YJKs colour-colour diagram 
and variability 

We showed above that most known QSOs lie in regions A+B 
of the YJK^ colour-colour diagram. However, the many other 
objects in the same region of colour space at first sight implies 
that spectroscopic follow up of the sample will produce only 
~ 6% which are actually QSOs, but we can also use the known 
A's variability over timescales of 40 - 600 days for 75% of the 
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200 400 600 



day 

Fig. 6. Ks band variation of a QSO with respect to its mean 
magnitude as indicated in Table |2] Nightly averaged points are 
connected by a continuous line while a dashed line represents 
a linear fit through them. Error bars are the standard deviations 
from the mean. The slope of the line is 1.8 +0.4x10"'* mag/day 
where day coiTesponds to the first observation in the VMC 
data. 



known QSOs to select a much tighter sample of candidates 
QSOs, a much greater fraction of which would be expected to 
be real QSOs when spectroscopically followed up. 

As 75% of the known QSOs show a variation in the 
band with a slope > 10""* mag/day (Fig.|7]) we expect 10 and 
23 QSOs with this characteristic in tile LMC 6_6 and 8_8, re- 
spectively, among the magnitude, error and quality flag lim- 
ited samples. These numbers have been confirmed by examin- 
ing the light-curve variation of sources in these two tiles, for 
which all Ks epochs have been obtained by the VMC survey 
over a time range of 300 - 600 days. 

To select objects with variability and colours like the 
known QSOs, we use criteria: photometric errors < 0.1 mag, 
quality flags = in YJKs, a slope of variation > 10 
mag/day, Y < 19.32 mag, J < 19.09 mag and K, < 18.04 
mag. We find a sample of 22 objects from tile LMC 8_8 of 
which 14, 7 and 1 located in region A, B and C respectively. 
This number is similar to the expected number of candidate 
QSOs derived earlier, i.e. 30 candidate QSOs are expected 
(Sect. [33T| l of which 75% (~ 23) show a variation in the K^^ 
band with a slope > 10"^ mag/day. The same criteria applied 
to tile LMC 6_6 result in a sample of 26 objects, 15 and 1 1 in 
regions A and B respectively, while the large number of ob- 
jects in region C reflects the contamination by reddened LMC 
variable stars. See Table H] for details. This suggests that the 
combination of YJK^ colours and variability selection is 
very effective in making a candidate QSO selection. The num- 
ber of candidate QSOs would be larger if all sources of photo- 
metric class A are used. 
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Fig. 7. Number of QSOs as a function of slope of variation in 
the band. Each bin is 0.5 x 10"^ mag/day. A vertical dashed 
line is drawn at 10""* mag/day; 75% of the QSOs have slope 
larger than this value. 

3.5.3. Identification of selected candidate QSOs 

A detailed discussion of the true nature of the individual QSO 
candidates is beyond the primary aim of this study and it 
is postponed to when spectroscopic data becomes available. 
QSO contaminants are discussed in detail in the next section. 
However, counterparts to the 22 and 26 candidate QSOs in 
tile LMC 8_8 and 6_6, respectively, have been searched for 
in SIMBAD. One match, an X-ray sources was found in tile 
LMC 8_8. In tile 6_6 there are 3 X-ray sources and 1 young 
stellar object (YSO) in region A, and 1 X-ray source and 2 
YSOs in region B. X-ray counterparts are usually found at a 
distance of up to 1 ' while YSO counterparts are found at dis- 
tances < 1". 

4. Discussion 

We have established YJK^; colour-colour and criteria to 
select candidate QSOs from VMC data. We now discuss the 
number and nature of objects which are not QSOs and which 
will also be selected by our criteria, to evaluate the degree of 
contamination of our sample of candidate QSOs that will not 
turn out to be QSOs. 

4.1. QSO candidate contaminants 

Photometric criteria for selecting QSOs often result in the mis- 
classification of sources that show similar colours and magni- 
tudes. The main contaminants in candidate QSO samples are 
normal galaxies, YSOs, PNe and evolved stars. Brown dwarfs 
have (Y - J) > 0. 8 mag while QSOs have (Y - J) < 0.8 
mag (Warren et al. I2007I I. The cool T-type dwarfs have also 
-1.5 < {J - Ks) < 0.5 mag (Burningham et al. |20T0] l while 
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L-type dwarfs have (7 - K^) > 0.5 mag (Knapp et al. l2004l l. 
Only the latter may be present in QSO regions BC (Fig. |2]l. 
In the next sections, the contribution by QSO candidate con- 
taminants is discussed in more detail using the photometric 
study by Gruendl & Chu ((2009 ), the spectroscopic studies by 
Kozlowski et al. (: 2012| l. Woods et al. (.201 I j and Kamath et 
al. (in prep.), and the recent study Kim et al. (I2012I I based on 
variability and photometric diagnostics. 

4.1.1. Gruendl & Chu ||2009) 

A statistically significant sample of QSO contaminants was 
assembled by Gruendl & Chu (2009) who analysed spectral 
energy distributions, from optical to mid-lR wavelengths, and 
mid-lR images from Spitzer in the central 7° x 7° area of the 
LMC. Their samples is: 949 definite and 126 probable back- 
ground galaxies, 855 YSOs, 56 PNe, 159 diffuse sources, 117 
evolved stars and 291 stars with IR excess have all been po- 
sitionally matched with the present VMC data. We have ex- 
cluded their lists of probable and possible YSOs since it is 
indicated that the first group may be background galaxies and 
the second group cannot be definitely classified. 

The distribution of the matched sources in the colour- 
colour diagram is shown in Fig. [8] There are: 286 definite and 
48 probable galaxies, 396 YSOs, 25 PNe, 123 diffuse sources, 
44 evolved stars and 165 stars. The majority of these sources 
have counterparts within ~ I" compared to an angular reso- 
lution of Spitzer of ~ 2" (Meixner et al. 2006 ). As expected, 
galaxies populate region A and B but most of them are located 
around the line separating them from the main locus of stars. 
This is also the preferred location of YSOs which extend to 
bluer colours and into region B. Most PNe are well confined 
in the region defined in Paper 11 while stars and diffuse sources 
occupy mostly regions C and D. A possible AGB star in the 
PNe region was found to be a PN by Miszalski et al. (201 la) 
and its colours are (F - 7) = 0.36 and (J - K^,) = 2.16 mag. 

The sample of background galaxies by Gruendl & Chu 
(120091 ) contains a mixture of marginally resolved extended 
sources and extended sources with a point source nuclear re- 
gion which were obtained from selection criteria tuned to find 
YSOs. They are therefore not representative of the background 
galaxy population in general. Among the sample of definite 
galaxies within the PN region Kamath et al. (in prep.) finds 
one QSO and one YSO out of three sources in common; the 
third one could not be classified due to low S/N spectral fea- 
tures. These findings suggest that the Gruendl & Chu's classi- 
fication, which is based on photometric criteria, may not be re- 
liable. Sources with a slope of variation in the band > 10 
mag/day represent a subset of the samples of QSO contam- 
inants and their distribution in the colour-colour diagram is 
similar to that of the full samples (Fig.|8]). 

By considering only sources with photometric errors < 0.1 
mag and quality flags = in each wave band, all YSOs are 
excluded (including the one above), only one star, a diffuse 
source and a PN remain in region B as well as two stars in re- 
gion C while the number of background galaxies is also con- 
siderably reduced and confined within regions A and B, but 
for one source (Fig.lS). The main criterion responsible for the 
removal of sources is the requirement on the quality flags. In 
fact, the magnitudes of most sources were extracted following 
a de -blending process corresponding to flag = 16 (the con- 
firmed QSO above has indeed this quality flag value in each 
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Fig. 8. Distribution of the Gruendl & Chu ( 120091 ) sources in 
the colour-colour diagram from present VMC data, (top row) 
All positionally matched sources, (middle row) Sources with 
a slope of variation in the band > lO""* mag/day. (bottom 
row) Sources with photometric errors < 0.1 mag and qual- 
ity flags = in the Y, J and Kf^ bands, (left column) Definite 
(blue circles) and probable (cyan triangles) background galax- 
ies, (middle column) YSOs (red circles), (right column) PNe 
(green squares), diffuse sources (yellow reversed triangles), 
evolved stars (magenta circles) and other stars (black trian- 
gles). Lines are as in Fig.|2l 



VMC wave band). Since the Gruendl & Chu's sample is un- 
confirmed, the remaining galaxies are either QSO contami- 
nants or QSO candidates themselves while the influence of 
other sources is negligible. These galaxies are not previously 
known, all but one have a well defined star-like appearance 
in the VMC images and about half of them have a slope of 
variation in the band > lO"'* mag/day. 

4.1 .2. Kozlowski et al. f20T2), Woods et al. i20m and 
Kamath et al. (in prep.) 

The spectroscopic investigation by Kozlowski et al. (120121 ) 
in the LMC produced not only a sample of 169 spectro- 
scopically confirmed QSOs but also a sample of non-QSOs 
comprising: 112 YSOs, 17 PNe, 39 Be stars, 24 blue stars 
and 66 red stars. They also observed that 390 sources with 
I > 19 mag are likely QSO candidates whose spectra are 
featureless and prevent a secure identification. Several non- 
QSOs were previously known. These objects were positionally 
matched to the VMC merged catalogue and counterparts were 
assigned to the nearest object within 1". This value is consis- 
tent with the mean distance to the QSOs derived in Sect. 12.31 
The Kozlowski et al.'s non-QSO sample was complemented 
with the Spitzer spectroscopically confirmed LMC sources by 
Woods et al. ( [20TlT ). These include: 90 AGB stars, 23 post- 
AGB stars, 7 PNe, 29 YSOs, 19 red supergiant stars and 7 
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background galaxies. Presently available VMC counterparts 
are at a distance of ~ 1". A further non-QSO sample con- 
sisting of red-shifted normal and Seyfert galaxies, spectro- 
scopically confirmed by Kamath et al. (in prep.) and cover- 
ing sources behind both the LMC and SMC, was also posi- 
tionally matched to the VMC data using the same criteria and 
included in the following discussion. The three studies consid- 
ered here produce a large sample of non-QSOs. This sample 
is large enough to permit a reliable estimate of the contamina- 
tion of QSO candidate samples by non-QSOs when the QSO 
candidates are selected on the basis of their position in the 
VMC colour-colour diagram and the slope of variation in the 
A's band. It is not the aim of this investigation to discuss the 
completeness of the spectroscopic non-QSO sample. 

The distribution of spectroscopically confirmed non-QSOs 
in the colour-colour diagram is shown in Fig.|9l There are 106 
galaxies of which 7 are Seyferts, 41 YSOs, 9 PNe, 14 post- 
AGB stars, 1 1 blue stars (including Be stars) and 37 red stars 
(including AGE and red super giant stars). Ninety-four likely 
QSOs are also matched with the VMC data and shown in Fig. 
|9] Comparing Fig.|2]and|9] it seems that QSO candidates se- 
lected from the left side of the green line in regions A and B 
are unlikely to be contaminated by stars, galaxies or YSOs, 
but there will be a contamination by PNe. QSO candidates se- 
lected from region A but to the right of the green line will be 
contaminated by galaxies and YSOs. QSO candidates selected 
from region B but to the right of the green line will be contami- 
nated by normal and Seyfert galaxies, YSOs and red stars. The 
sub-set of variable objects (with slope in the band > 10""^ 
mag/day) weeds out contaminating red stars from regions B 
and C, and PNe from regions A and B while the distribution 
of galaxies and YSOs is unchanged. The distribution of ob- 
jects with photometric errors < 0.1 mag and quality flags = 
results in a negligible number of galaxies and YSOs in regions 
B and C while none remain in region A. The Kozlowski et al.'s 
candidate QSOs populate the colour-colour diagram similarly 
to known QSOs (Fig. |2]i and those that remain after applying 
variability, photometric uncertainty and quality flag criteria, 
are high confidence candidates also according to our study. 

4.1.3. Kim et al. (20m 

The MACHO database has been thoroughly mined for the 
identification of high confidence QSO candidates behind the 
LMC by Kim et al. (2012). Initially, using a method purely 
based on the QSO time variability, 2566 candidates were iden- 
tified (Kim et al. |201 11 1. With the additional information from 
observations at other wavelengths and the development of a 
model based on the diagnostics of known QSOs, 663 high con- 
fidence QSOs were later identified (Kim et al. |2012| l. 

Figure [To] shows the distribution of the Kim et al.'s can- 
didate samples positionally matched with the VMC data in 
the colour-colour diagram. There are 574 low confidence and 
179 high confidence QSOs with a VMC counterpart within 
1 .5" in the merged catalogue. The major point to notice is that 
low confidence QSOs are distributed mainly in region D while 
high confidence QSOs occupy mostly regions A and B. This 
is the case in general, for sources with a slope of variation in 
the band > 10"^ mag/day and for sources with photomet- 
ric errors < 0.1 mag and quality flags = in each wave band 
strongly supporting the criteria developed in this study. 
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Fig. 9. Distr ibution of the Kozlowski et al. ( |2012| l. Woods 
et al. ( 1201 11 1 and Kamath et al. (in prep.) spectroscopically 
observed sources in the colour-colour diagram from present 
VMC data, (top row) All positionally matched sources, (mid- 
dle row) Sources with a slope of variation in the band 
> 10"** mag/day. (bottom row) Sources with photometric er- 
rors < 0.1 mag and quality flags = in the Y, J and bands, 
(left column) Likely QSOs (black circles), normal (blue cir- 
cles) and Seyfert galaxies (cyan triangles), (middle column) 
YSOs (red circles), (right column) PNe and post-AGB stars 
(green squares), red stars (magenta circles) and blue stars 
(black triangles). The lines are as in Fig.|2] 



4.2. QSO variability 

While the optical variability of QSOs is well established (e.g. 
MacLeod et al. 20 11) much less is known in the near-IR. The 
most recent investigation by Kouzuma & Yamaoka ( 120121 ) ad- 
dresses the properties of the ensemble variability of QSOs 
using two data points from the Two Micron All Sky Survey 
(2MASS; Skrutskie et al. |2006l l and the DEep Near-Infrared 
Southern sky survey (DENIS; Epchtein et al. 1999 ) or the 
UKIR T Infrared Deep Sky Survey (UKIDSS; Lawrence et al. 
|2007] l. About 15 - 25% of the QSO emission in the K, band is 
attributed to the accretion disc (Tomita et al. 2006) and the re- 
maining to thermal radiation by dust in the AGN torus, heated 
by the disc (Barvainis |1987l l. The reason of variations in the Ks 
emission may be due to changes in the inner radius of the torus 
as a consequence of variations in the UV/optical flux from the 
AGN (Koshida et al. 2009 Kawaguchi & Mori 201 1 ). 

The best sampled near-IR QSO light-curves were ob- 
tained by the Multicolour Active Galactic Nuclei Monitoring 
(MAGNUM) programme on NGC 4151 and NGC 5548 
(Yoshii, Kobayashi & Minezaki |2004l Kosh ida et al. |2009] l or 
focused on 3C273 (Courvousier et al. 1 19881 ) for a total of three 
sources. In the first study the K light-curves of the AGN nuclei, 
after subtracting the contribution from the host galaxy, show 
clear minima and maxima each over a time range of ~ 200 



14 



Cioni et al.: QSOs behind the Magellanic system 



2 
1.5 

0.5 

2 

1.5 
^" 1 

0.5 

2 

1.5 

f' ^ 
0.5 





•A V ' 

.■ D 



0.5 

Y-J 



- D ^ 

1 1 1 1 1 1 1 1 1 1 1 1 1 ~ 


A P 

"^-tx • ■ 

D 

' 1 1 1 1 1 . 


1 1 1 1 J 


; ? 

D 

' 


• c - 



0.5 

Y-J 



Fig. 10. Distribution of the Kim et al. ( |2012| l sources in the 
colour-colour diagram from present VMC data, (top row) All 
positionally matched sources, (middle row) Sources with a 
slope of variation in the band > 10"** mag/day. (bottom 
row) Sources with photometric errors < 0.1 mag and quality 
flags = in the Y, J and K^^ bands, (left column) Low con- 
fidence QSOs (blue circles), (right column) High confidence 
QSOs (magenta circles). Lines are as in Fig.|2] 



days. Overall the K variation is smoother than the variation in 
the V band which appears irregular and associated with rapid 
variations on time scales of several days. The luminous 3C273 
source shows fast and structured flaring events also in the near- 
IR and only ~ 15 days apart. Near-IR variations of QSO emis- 
sion over a time range of up to twenty years but with a sparse 
sampling were also found by Neugebauer et al. (11989l l. 

The VMC multi-epoch observations represent a consider- 
able step forward in the study of QSO near-lR variations. The 
present sample shows both smooth trends, i.e. an increasing 
and decreasing brightness over a time range of 300 - 600 days, 
and structured events, i.e. broad (up to 100 days wide) and 
narrow (20 - 40 days wide) peaks, in agreement with previous 
studies. Despite the large number of QSO light-curves anal- 
ysed in this study, there is no clear pattern of variation, but it is 
possible that similarities in the QSO light-curves occur across 
a time range of several years. In the future, when more QSO 
light-curves are available across a wide time range, we plan 
to construct empirical template light-curves from the combi- 
nation of individual QSO light-curves with similar shape and 
structure functions (e.g. Hughes, Aller & Aller ,1992) . These 
can then be used to refine the search for QSOs and eventually 
to derive their redshift as in Dai et al. (2012). A theoretical 
interpretation of the VMC near-lR light-curves and structure 
functions may also provide further evidence for the mecha- 
nism responsible for the variations and contribute to studies of 
the structure and dynamics of the AGN (e.g. Wold, Brotherton 
& Shang 120071 Ovcharov et al. |2008] l. 



4.3. QSO numbers behind the Magellanic Clouds 

Assuming that QSOs are homogeneously distributed across 
the surface of the Magellanic system we can use the number of 
VMC candidate QSOs derived in the previous sections to esti- 
mate the total number of QSOs detectable by the VMC survey. 
Using the magnitude limited sample in tile LMC 8_8, which is 
not aff'ected by interstellar reddening and crowding since it is 
located in the outer disc of the galaxy, we derived that 22 ob- 
jects are candidate QSOs. The VMC survey comprises of 1 10 
tiles of which 68 in the LMC, 27 in the SMC, 13 in the Bridge 
and 2 in the Stream areas. The expected QSO population ac- 
cessible to the VMC survey is thus of at least 1500 objects 
behind the LMC, 600 behind the SMC, 300 behind the Bridge 
and 50 behind the Stream areas. This means that ~ 85% of 
the VMC accessible QSO population behind the Magellanic 
system is yet to be found. In addition, these numbers are con- 
servative limits since they have originated from the application 
of magnitude, error, quality flag and variability restrictions. 



5. Conclusions 

The availability of large samples of QSOs across the 
Magellanic system is of fundamental importance as an astro- 
metric reference source for investigations of the proper motion 
and for a study of the interstellar and intergalactic medium 
along the line of sight as well as to address the nature of the 
extra-galactic sources themselves. The search and characteri- 
sation of QSOs behind the Magellanic system is an on-going 
process as part of follow-up observations of large samples se- 
lected from optical light-curves and multi-wavelength pho- 
tometry (Kozlowski et al. 2012) . For a most efficient use of 
telescope time the determination of high confidence QSO sam- 
ples is necessary (Kim et al. l2012l l. 

We present criteria for selecting candidate QSOs based 
solely on the multi-epoch near-lR photometry from the VMC 
survey developed from the properties of 1 17 known QSOs ob- 
served behind the Magellanic system. Most known QSOs oc- 
cupy two neighbouring regions of the (Y - J) vs. (J - Ks) di- 
agram, one is dominated by QSOs with a star-like appearance 
and the other by QSOs with a galaxy-like appearance. The 
analysis of multi-epochs shows that ~ 75% of the QSOs 
have a slope of variation > 10 mag/day across a time range 
up to 600 days. The combination of YJK^ colour-colour and 

variability criteria, with appropriate data quality flags, can 
be used to define samples of candidate QSOs which should be 
relatively free of contamination by non-QSOs. 

QSO candidates in LMC tile 6_6, including the 30 Doradus 
region, and 8_8, including the South Ecliptic Pole region, 
have been extracted from sources with photometric errors 
< 0.1 mag and quality flags = in each VMC wave band 
and brighter than 19.32, 19.09 and 18.04 mag in the Y, J 
and band, respectively. These magnitudes correspond to 
the faintest QSO in the sample by Kozlowski et al. ( I2012I I. 
Thirteen and 22 objects, respectively, satisfy the colour-colour 
and variability criteria. The selection of QSOs from the VMC 
data is in excellent agreement with the sample of high confi- 
dence QSOs by Kim et al. (|20T2] |. 

The contamination by non-QSOs has been estimated us- 
ing the photometric samples of YSOs, PNe, diffuse sources, 
evolved stars and other stars with IR excess, from Gruendl & 
Chu ( 120091 1 and spectroscopic samples from Kozlowski et al. 
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(I2OT2I 1. Woods et al. dSOTlT l and Kamath et al. (iti prep.). PNe 
represent the major source of contamination in the regions of 
the colour-colour diagram dominated by star-like QSOs while 
normal galaxies and YSOs are the major contaminants in the 
region dominated by galaxy-like QSOs. The sub-set of sources 
with a slope of variation in the band > lO""* mag/day is not 
influenced by PNe and red stars. The sub-set of sources with 
photometric errors < 0.1 mag and quality flags = in each 
VMC wave band is negligible. Since known QSOs may be 
de-blended and have quality flags =16 the latter represents a 
reliable whilst not complete sample of candidates. 

VMC magnitudes and light-curves for known QSOs 
are provided in a table and in the appendix. The full sample 
of known QSOs comprises 332 QSOs from the literature and 
28 newly confirmed QSOs by Kamath et al. (in prep.). VMC 
data on a fully observed tile in the outer LMC disc (8_8) have 
been used to estimate the number of QSO candidates behind 
the entire system as covered by the VMC survey. This is of 
the order of 1500 behind the LMC, 600 behind the SMC, 300 
behind the Bridge and 50 behind the Stream areas. Compared 
to the known number of QSOs in these regions ~ 85% of the 
population is still missing. The VMC survey is the most sen- 
sitive near-IR survey of the Magellanic system to date pro- 
viding for the first time near-IR counterparts to many QSOs. 
Spectroscopic observations of the candidates identified here 
will test the method and extend the selection to faint magni- 
tudes exploiting the whole VMC range available. 
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Appendix A: band light-curves of known QSOs 

Figure I A. II shows the band light-curves of known QSOs 
in the presently available VMC data (VSA data release 
VMCv20120126) where single epoch data are included for 
completeness. 
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Fig. A.l. Ks band light-curves of known QSOs with respect to their mean magnitude. From top to bottom and left to right they 
follow the order in Table |2] Points represent the simple average of observations obtained within the same night/epoch and error 
bars correspond to the error of the mean. A continuous line connects the data points while a dashed line represents a linear fit 
through them. Day corresponds to the first observation in the VMC data. Single epoch data are included for completeness. 
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Fig. A.l. continued. 
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Fig. A.l. continued. 
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Fig. A.l. continued. 
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Fig. A.l. continued. 
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